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Abstract 

We formulate the framework of M- fold supersymmetry in one-body quantum mechanical systems 
with position-dependent mass (PDM). We show that some of the significant properties in the 
constant-mass case such as the equivalence to weak quasi-solvability also hold in the PDM case. 
We develop a systematic algorithm for constructing an A^-fold supersymmetric PDM system. We 
, apply it to obtain type A AA-fold supersymmetry in the case of PDM, which is characterized by 

Q \ the so-called type A monomial space. The complete classification and general form of effective 
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potentials for type A AA-fold supersymmetry in the PDM case are given. 
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I. INTRODUCTION 



Recently, quantum systems with position-dependent mass (PDM) have attracted much 
attention in various research fields of physics such as semiconductors, quantum dots, liq- 
uid crystals, and so on. Accordingly, investigations into exact solutions of PDM quantum 
systems have been carried out increasingly in the last few years. For references, see e. 



Refs. QiliaBQiailinBHQHHH 

and those cited therein. Due to the fact that the position-dependent mass m(q) does not 
commute with the momentum operator p = — id/dg, ambiguity arises in defining a quan- 
tum kinetic operator which is formally Hermitian and reduces to the classical kinetic term 
T = p 2 /2m(q). Hence, the following operator proposed by von Roos has been generally 
considered: 



H = — ( m{q) a ■^-m{q) 13 ■^-m{q) 1 + m(g) 7 -^-m(g) /3 — m{q) a 
4 \ dq dq dq dq 



V(q), 



with the constraint a + /3 + 7 = —1. A different choice of the parameters results in a different 
correction to the original potential profile V(q), and the above Hamiltonian always has the 
following form: 



H 



2m(q) dq 2 2m(q) 2 dq 
where the effective potential U (q) is given by 

tt, \ it/ \ , \ m"(q) / \ m'(q) 2 

U(q) = V(q) - (a + 7) - + (try + a + 7)- 



;i.2) 



4m(g) s 



2m(q)' c 



;i.3) 



Hence, the typical investigations into (quasi-)exact solvability of PDM quantum systems 
consist in finding simultaneously a pair of an effective potential U(q) and a mass function 
m(q) for which the PDM Hamiltonian (jl.2j) admits (a number of) exact eigenf unctions in 
closed form. Up to now, two different methods have been frequently employed, namely, 
coordinate transformations including point canonical transformations 1 11 |2j. la l a 171,181 l9l ll2L 
0E10@, and supersymmetric methods 0, B 0, H Q E3, 13 EJHI2JI2I Ul Ef 
The latter approaches were also applied to many-body PDM quantum systems |26j | . For the 
methods themselves developed in ordinary quantum systems, see references cited therein. 

On the other hand, the framework of A/"-fold supersymmetry has been rapidly devel- 
oped in one-body ordinary quantum mechanical systems, especially in the last few years. 
It was originally proposed as a higher derivative generalization of ordinary supercharges 
(28l | . Later, a significant breakthrough was achieved by the proof of the equivalence between 
A/"-fold supersymmetry and (weak) quasi-solvability [29]. Based on the equivalence, a sys- 
tematic algorithm for constructing an A/"-fold supersymmetric system was developed [30]. 
Owing to the fact that the framework of A/"-fold supersymmetry includes all the ordinary 
supersymmetric methods as its special cases and to its equivalence to weak quasi-solvability, 



which is less restrictive concept than quasi-exact solvability 31, 32], it provides one of the 



most powerful methods for investigating one-body quantum mechanical systems which admit 
analytic solutions in so far the least restrictive sense. 

Considering the present situation described above, it is natural to ask whether or not 
we can formulate A/"-fold supersymmetry also in PDM quantum systems, especially as a 
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powerful tool for constructing (quasi-) solvable PDM systems. In this paper, we show that 
it is indeed possible and that the characteristic properties such as the equivalence to weak 
quasi-solvability also hold in the PDM case. Furthermore, we generalize the systematic 
algorithm for constructing an TV-fold supersymmetric ordinary quantum system in Ref. [3G| 
to the PDM case. As an illustration, we apply it to construct so called type A A/"-fold 
supersymmetry j3j| |34| in PDM quantum systems. 

The paper is organized as follows. In the following section, we review the precise definition 
of some important concepts in the paper, such as quasi-solvability, to avoid ambiguity. In 
Section IIII1 we define A/"-fold supersymmetry in one-body quantum systems with position- 
dependent mass and discuss a couple of general significant consequences. In Section IIV1 we 
develop a systematic algorithm for constructing an Af-iold supersymmetric PDM quantum 
system by slightly generalizing the one in the constant-mass case in Ref. |3Cj. In Section IVl 
we apply the algorithm to construct type A A/"-fold supersymmetry in the PDM case. In 
Section EH we completely classify and present the explicit forms of all the inequivalent type 
A models with arbitrary position dependence of mass. In the last section, we summarize 
the results and discuss various future issues. 



II. DEFINITION 

First of all, we shall give precise definition of some key concepts in the paper. The original 
ideas can be found in Refs. [ijj |35| but we have slightly modified the terms from the view 
point of the recent advances in this research field. Let if be a linear differential operator of 
a single variable and Pa/ be an (at most) A/"th-order linear differential operator. Then, H is 
said to be weakly quasi- solvable with respect to Pa/ jlH if it preserves the vector space V^r 
defined by 

V A r = kerP Ar . (2.1) 

A linear differential operator if of a single variable is said to be quasi- solvable if it preserves 
a finite-dimensional functional space W whose basis admits an analytic expression in closed 
form: 

HVxCVx, dimV M = n(N) < oo, V M = <0i(g), . . . , <j> n {N){<l)) ■ (2.2) 

It is evident that a weakly quasi-solvable operator is also quasi-solvable if linearly indepen- 
dent Af solutions of Pa/0 — can be obtained in closed form. An immediate consequence 
of the above definition of quasi-solvability is that, since we can calculate finite-dimensional 
matrix elements H^i defined by 

n(Af) 

H<p k =J2UkA, k = l,...,n{N), (2.3) 
i=i 

we can diagonalize the operator H and obtain its spectra in the space VV with finite algebraic 
manipulations. However, these calculable spectra and the corresponding vectors of W in 
general only give local solutions of the characteristic equation. This fact naturally leads 
to the well-known concept of quasi-exact solvability. A quasi-solvable operator H is said 
to be quasi- exactly solvable (on S C R or C) if the invariant space W is a subspace of a 
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Hilbert space L 2 (S) on which the operator H is naturally defined. It is evident that if an 
operator is not only quasi-solvable but also quasi-exactly solvable, the calculable spectra 
and the corresponding vectors of W gi ye a P ar t of the exact eigenvalues and eigenfunctions 
of H, respectively. On the other hand, there are cases in which all the eigenfunctions and 
eigenvalues can be obtained analytically such as the well-known harmonic oscillator. In order 
to characterize these cases, we first introduce another subclass of quasi-solvability. A quasi- 
solvable operator H is said to be solvable if it preserves an infinite flag of finite-dimensional 
functional spaces W, 

Vi C V 2 C ■ ■ • C V N C • • • , (2.4) 

whose bases admit explicit analytic expressions in closed form, that is, 

HV N dV N , timVx = n{N) < oo, = (^(g), . . . , n (Ao(?)>, (2-5) 

for all M = 1, 2, 3, . . .. A consequence of solvability defined above is that, for an arbitrary 
natural number A/", we can obtain additional n(N '+ 1) — n(N) local solutions of the charac- 
teristic equation for H in VV+i with finite algebraic manipulations, based on the knowledge 
of them in the subspace Vjv C VV+i- Hence, arbitrary number of local solutions of the char- 
acteristic equation for H can be calculated, in principle, with a well-defined finite algebraic 
algorithm. As in the case of quasi-solvability, however, they are not necessarily the exact 
eigenfunctions and eigenvalues of H . A solvable operator H is said to be exactly solvable 
(on S C K or C) if the sequence of the spaces (12. 4 j) in S satisfies, 



V N {S)^L 2 (S) (Af^oo). (2.6) 

This definition precisely fits what is commonly meant by exactly solvable. In the PDM case, 
it is pointed out ji^] that in addition to the square integrability, any eigenfunction ip of a 
PDM Hamiltonian (jl.2|) should satisfy on the boundary dS 

^ - (, ^ 8S), (2.7) 
m(q) 2 

in order for the PDM Hamiltonian to be Hermitian. In this paper, however, we do not discuss 
Hilbert spaces, boundary conditions, normalizability, etc. which strongly rely on a specific 
choice of m(q). We do not assume Hermiticity of Hamiltonians and thus reality of effective 
potential U(q) and mass function m(q) either, taking into account possible application to 
recently developing non-Hermitian quantum theories. 

III. AA-FOLD SUPERSYMMETRY IN PDM QUANTUM SYSTEMS 

To define A/"-fold supersymmetry in one-body PDM quantum systems, let us first intro- 
duce fermionic variables ip and ip^ satisfying 

{^,^} = {^,^} = 0, {^,^} = 1. (3.1) 
In this paper, we consider a super-Hamiltonian 

T1^ = H^ + H +^ (3 .2) 
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where the components Hj^ are given by PDM Schrodinger operators: 



Hm ~ 2m( q )dq* + 2m(qydq + UAq) - ^ 
For the above system, we define TV-fold supercharges Q^- by 

Qy = Pjftf = Qm = Pfrl> = Pfrl>, (3.4) 

where the operator Pjy is an A/th-order linear differential operator of the following form: 1 



In Eq. (|3.4|) . the superscript t denotes the formal transposition. A system EL\/ is said to be 
M-fold supersymmetric with respect to if it commutes with them: 

[Q±,EEy]=0. (3.6) 

From the above simple generalization of A/"-fold supersymmetry to PDM systems, we can 
show that most of the relevant consequences of A/"-fold supersymmetry in ordinary quantum 
systems with constant mass also hold in the present case. Let us first prove the equivalence 
between A/"-fold supersymmetry and weak quasi-solvability. Suppose a super-Hamiltonian 
ELy is AA-fold supersymmetric. It is easy to see that the component operators Hj^ satisfy 
the following intertwining relation 

P M H M = P+P^, (3.7) 

and its formal transposition, P^-H^- = H^Pj^, and thus they preserve the vector spaces 
ker P^, respectively. Hence, are weakly quasi-solvable. Conversely, let us assume that a 
PDM Schrodinger operator H is weakly quasi-solvable with respect to an A/th-order linear 
differential operator Pa/ of the form (|3.5|) . Define a linear operator G = P^H — H\P^ with 

rr rr ^ , ^ ~ 1 ^ , ^^m!' 3A" 2 (m') 2 

From the assumed weak quasi-solvability, H ker Py C ker Pa/, we have 

0G ker P N G<f> = Pj^Hcj) - H x P U (j) = 0. (3.9) 

On the other hand, it is easy to check that G is of order (at most) A/" — 1 while the dimension 
of ker Pa/ is N . Hence, Eq. (13. 9J) cannot be hold unless G is a null operator. Therefore, if 
we put Hjj = H, Hjj = Hi, P^ = P//, and P^ = Pfc, they satisfy the intertwining relation 
(|3.7|) and thus compose an A/"-fold supersymmetric system. 

Another significant consequence of AAfold supersymmetry in ordinary quantum systems 
is that the anti-commutator of A/"-fold supercharges Qj^ and Q^- is a polynomial of degree 
M in the super-Hamiltonian ELy, the polynomial being (proportional to) the characteristic 



Here we fix the irrelevant overall multiplicative constant factor so that Py is monic when m(q) = 1. 
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poly nomial of the component Hamiltonians Hjj- restricted in the invariant subspaces Vjy 



{Q^, Q+ } cc det (F±| v ± - Ha,) . (3.10) 

We can easily prove that the same is also true in the PDM case. Indeed, in the proof of 
the theorem on SUSY algebras with T symmetry in Ref. [36j], we only need to modify the 
functions wj^- in Eq. (6) as wj^(x) = (=)=l) A/ 'm(x) _ ~ and the definition of the operators rf 
in Eq. (19) as 



r t =m(x) a + xi (x)J , = (^--^- + xt(x^Jm(x) a, l = l,...,Af, (3.11) 

for a proof of Eq. ()3.10|) in the PDM case. In our definition of Hamiltonians ()3.3j) . the 
proportional constant in Eq. ()3.10|) is 2^. 

Other characteristic features such as the almost isospectral property between a pair of 
Hamiltonians Hj^ are easily derived in the same way as in the ordinary case (cf. Ref. [29]), 
and thus we don't repeat a presentation of them in this paper. 



IV. ALGORITHM FOR CONSTRUCTING AN AA-FOLD SUSY PDM SYSTEM 

In Ref. |3(j , a systematic algorithm for constructing an A/"- fold supersymmetric ordinary 
quantum mechanical system was developed. In this section, we show that it can be easily 
generalized to the PDM case. The starting point is an A/"-dimensional linear functional space 

Vx={y x {z),...,m{z)), (4.1) 

and a second-order linear differential operator 

Hi=-A{z)^-B(z)±-C{z), (4.2) 
which leaves Vjv invariant. Let 




denote the most general AAh-order linear differential operator with kernel VV, where the 
function g(z) is for the time being undetermined. We shall first construct another second- 
order linear differential operator of the form 

H+ = H^-6C(z), (4.4) 

satisfying the intertwining relation 

P N H N - H^P N = 0. (4.5) 
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To this end, note that the l.h.s. of Eq. ()4.5|) is in general a linear differential operator of 
order M + 1. Equating to zero the coefficients of d^ +l and d?f in this operator, we obtain 
the following two equations for the functions g{z) and 5C(z): 

sc = m^. (,» -^y M (B'-™)- - 2 ^_, (4.7) 

When Eqs. (|4.6j) and (J4.7j) are satisfied, the l.h.s. of Eq. (J4.5|) is a linear differential operator 
of order at most A/"— 1 annihilating the M- dimensional vector space Vjj, and hence it vanishes 
identically. 

The last step in our construction consists in applying a change of variable 

z = z{q) (4.8) 

and a gauge transformation 

H±^e-^H±e w «\ z=z{q) = H±, (4.9) 

to simultaneously convert H^j- into the PDM Schrodinger form (J3.3)) . Note that it is certainly 
possible since (by construction) HJ^ and Hjj differ by a scalar function only. The appropriate 
change of variable and gauge transformation are determined by 

z\qf = 2m{q)A{z), (4.10) 
dW^ _ z"(q) m(q)B(z) m'{q) 



The effective potentials Ufr are given by 



Ut(q) 



dq 2z'(q) z'(q) 2m(q) 

JV 

1 



(4.11) 



2m(q) 



dW^V d 2 W^ , m\q)dW N 



dq J dq 2 ^ m(q) dq 



-C+(z(q)), (4.12) 



where C (z) = C(z) and C + (z) = C(z)+5C(z). From the above construction it immediately 
follows that the system (J3~2j) and (jS3j), with Hjfr given by Eq. (jOj) and P N by 

Px = e-^P^\ z=z{q) , (4.13) 

is TV-fold supersymmetric. Indeed, intertwining relation ()3.7|) follows by applying the gauge 
transformation and change of variable to Eq. (|4.5|) . Furthermore, it is important to note 
that the form of Pjj in Eq. ()4.13|) is compatible with Eq. (|3.5jl . From Eqs. ()4.6|) and ()4.10|) 

we have 

9'{z) = Afz"(q) _ Mm\q) 

g{z) z'{q) 2 2m{q)z'{q)' 1 ' ' 

Integrating the latter equation we obtain 

g(z) = m(q)-%z'(q) Af , (4.15) 



7 



where we take the proportional constant to be 1. Thus, it immediately follows from Eqs. ()4.3|) 
and ()4.13|) that P^ in Eq. (|4.13|) is indeed of the form ()3.5|) . It is evident from the construc- 
tion that the Hamiltonian Hjj preserves the kernel of Pjj given by 

V M = kei P M = (^(q), . . . , m (q)), (4.16) 

with 



<Pi(q) = e w x<pi{z) 



z=z(q) 



,AT). 



(4.17) 



We thus call the space V^- the solvable sector of Hjj. 

Although the construction of an A/"-fold supersymmetric PDM system itself has been 
completed, we can make it entirely symmetric with respect to the partner Hamiltonians H7f 
and H^f. For this purpose, note that from the transposition of the intertwining relation 
()3.7|) . it follows that leaves invariant the kernel of the supercharge 



(4.18) 



Using the identity (d z Y = —z'(q)d z z'(q) 1 and Eq. ()4.3|) with (|4.15|) . we can express P£ as 

P+ = e- w "P£e w ", (4.19) 

where 



P+ = m (q) 2 Z '( q ) 
and the function VV^ is given by 



d 

dzj 



fc=0 



d Y 



dzj 



(4.20) 



W7r+{N- l)ln|z'(g)| - — ln|m(g)| 



(4.21) 



Now, the partner Hamiltonians H^- are expressed in a completely symmetric way as 



H 



J = V U a7/ a .<. U A. 



(4.22) 



charges = kerP^ and V^- = kerP^, respectively. To express W^- in a symmetric way, 



where the gauged Hamiltonians and leave the kernels of the gauged AA-fold super 

= 1 

we introduce two functions by 

dtt# 



dq 



E{q) 



z"(q) 
z'(q) 



(4.23) 
(4.24) 



From Eq. (|4.10|) . its immediate consequence 

z"(q) = m(q)A'(z) + 



m'(q)A(z) 
z'(q) 



(4.25) 
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and Eq. ([4.110 . the function W(q) is expressed as 

m(q) ( 'N -2 



W(q) 



z'{q) 



A'(z) + B(z] 



m(q) 



Q(z). 



(4.26) 



We then have 



■A/", , / M Af ~ 1 
-— ln|m(g)| + 



!, , / m 
-\n\m(q)\ + 



dqE(q) T / dqW(q) 

m(q)Q(z) 



]n\2A(z)\ ± / dz 



2A(z) 



(4.27) 



The connection between the gauged Hamiltonians Htf and follows easily from Eqs. (J4.9[) . 
(JOS), and dOS|) as 



= exp 



dgW(g) if^exp 2 / dqW(q) 



(4.28) 



Using Eqs. (fOjl . ([Ojl . (I4~7I) . (j4~TU[l . and (j4~2E|) . we obtain the following unified formula for 
the gauged Hamiltonians: 



'Af - 2 



dz 



-C(*) 



-(1±1) 



A/"- 1 



Q'( z ) ~ ^A\z)wm-i(z) - A(z)w' M _ 1 (z) 



(4.29) 



Interestingly, the form of the gauged Hamiltonian Hj^ as well as Hj^- is completely the same 
as in the case of constant mass, cf. Eq. (2.45) in Ref. 113]. It should also be noted that 
the mass dependence of the Hamiltonians H^- given by Eq. (|4.22[) emerges only through the 
change of variable determined by Eq. (|4.1U[) and the gauge potentials by Eqs. (|4.11j) and 
(I4.21[) since the gauged Hamiltonians (j4.29[) do not depend on the mass function m(q). As 
a consequence, each of the spectrum of the Hamiltonians Hj^- does not depend on the mass 
function m{q)? This kind of spectral independence from the mass function is first pointed 
out in Ref. Efl. 



V. TYPE A AA-FOLD SUPERSYMMETRY 

We shall now construct type A A/"-fold supersymmetric PDM quantum systems with the 
aid of the algorithm just developed in the previous section. Type A M-io\d supersymmetry 
|33l ] is characterized by the so-called type A monomial space (3Gj preserved by Hj^: 



Vm=(1,z,...,z^- 1 ). 



(5.1) 



Applying the algorithm with this type A monomial space, we can construct the most general 
form of type A A/"-fold supersymmetric PDM quantum systems. Fortunately, we can omit 



2 Some of the eigenvalues can (dis)appear depending on m(q) since normalizability of the corresponding 
wave functions does rely on it. 
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the process of constructing gauged Hamiltonians. Noting the fact that the form of gauged 
Hamiltonians iJjy is completely the same as in the constant-mass case, we immediately have 



d 2 



J ^-^A'(z)±Q(z) 



d 

dz 



12 2 
where R is a constant, and Q(z) and A(z) must satisfy 

d 3 g(^ x 



(5.2) 



dz 3 
d 5 A(z) 

dz 5 



for Af > 2, 
for TV > 3, 



or equivalently, 



Q(z) = b 2 z 2 + b x z + b 



for Af > 2, 



= a 4 z 4 + a 3 z 3 + a 2 z 2 + a\z + a for A/" > 3, 



(5.3) 
(5.4) 



(5.5) 
(5.6) 



where hi and a« are constants. The most general A/th-order linear differential operator of the 
form ()4.3|) . with g(z) given by Eq. f)4.15|) . which annihilates the type A space is obviously 



.V 



P N = m (q) 2 z '( q ) 



d" 
dz*' 



(5.7) 



Hence, using Eqs. fl4.13|) . ()4.24j) . and ()4.27j) . we obtain the operator P^ for the type A 
A/"-fold supercharge 



P M = m (q)-% Yl ( 

fc=0 \ 



d Txr . . Afm'(q) Af - 1 - 2k . 

— + W(q) - - ;\ J + E{q) 

dq Am(q) 2 



(5.8) 



where the products of operators are ordered according to the following definition: 

fci 



k=ko 



A 



fco • 



(5.9) 



The above Pjj indeed reduces to the ordinary type A AA-fold supercharge when m(q) = 1. 
The pair of type A PDM Hamiltonians Hj^- can be expressed in terms of the functions E(q), 
W(q), and m(q) using the following formulas, which follow from Eqs. ()4.10j) and ()4.26|) : 



A'(z) 
A"(z) 
Q'(z) 



z(q) 



fE{q) m'(q) 
\m(q) 2m(q) 2 
E'(q) + E{qf 3m'{q)E{q) m{q)m"(q) - 2m'{qf 
m(q) 2m(q) 2 2m(q) 3 

W'(q) + E(q)W(q) m'(q)W(q) 



m(q) 



m(qy 



(5.10) 
(5.11) 
(5.12) 
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From Eqs. ()4.22j) . (|4.27jh (|5.2|h and the above formulas, we finally obtain 

H ± = _ 1 d 2 + m^)_ d_ + W[qf _ J\f 2 - 1 , _ 2) 
^ 2m(q) dq 2 2m{q) 2 dq 2m(q) 24m(g) 

A/~ 2 + 2m"(g) _ 5Af 2 + 16 m'(q) 2 f W'jq) _ m'{q)W{q) \ 

24 m(g) 2 96 m(g) 3 VM?) 4 ™(<?) 2 / 

This is the most general form of type A A/"-fold supersymmetric PDM Hamiltonians and 
exactly reduces to the ordinary type A when m(q) = 1. It should be also noted that the 
above formula is consistent with Eq. ()3.8|) . From the form of the type A M- fold supercharge 
f!5.8|) we have 

M Af 2 JV+2 

Wj\f-i(q) = Mm(q)~^W(q) — m{q) ~m!(q). (5.14) 

With this formula, it is straightforward to check for the type A Hamiltonians that 
+ _ _ AfW _ Afm'W 

u-2 , M-l m-a , M 2 m" 3Af 2 m' 2 



m 2 w A f_ 1 -\ — m 2 mw A /-_ 1 + 



4m 2 8m 3 



and thus consistent with Eq. ()3.8|) . Finally, we shall express the two conditions ()5.3j) and 
(15. 4 j) for type A M-fo\d supersymmetry in terms of E(q), W(q), and m(q). Using Eqs. ()4.1U|) . 
KZQ . and (jOBj) . we obtain 

±-e)±(±+e)™ = for AT>2, (5.15) 
dq J dq \dq J m 

±-2E)(± -e] A (± + E ) (B_ 4),o for Af>3. (5.16) 
dq J \dq J dq \dq J \m 2m 2 / 

In particular, in the case of M = 1 the system composed of Pa/ and reduces to 

m (q)h \dq 4m(q)J 
, 1 d 2 m'(q) d W(q) 2 m"(q) 7m'(q) 2 



2m(q) dq 2 2m(q) 2 dq 2m(q) 8m(q) 2 32m(g) 3 

± M_«))_ B , ( , 18 ) 

with which the super-Hamiltonian H x and supercharges defined by Eqs. (|3.2|) and (|3.4|) 
satisfy ordinary superalgebra: 

[Qf,Hj=0, {Qr,Ql} =2(H M + R). (5.19) 

Hence, it exactly reduces to (ordinary) supersymmetric PDM quantum systems. We note 
that thanks to the assignment of the functions W(q) and m(q) in Eq. (|5.17jl . which is 
different from the one conventionally employed in the literature, the resulting pair of H l 
has completely the symmetric form. 
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Case 


Canonical Form 


f(u) 


T 

1 


1/2 


u 


II 


2z 


u 2 


III 


2vz 2 




IV 


2v(z 2 - 1) 


cosh 2^[vu 


V 


2z 3 - g 2 z/2 - g 3 /2 





TABLE I: Canonical forms of A{z) and the functions f(u) which characterize the change of variable. 
The parameters v, g 2 , g% G C satisfy v / and g\ — 27 g 2 / 0. 



From Eqs. ()4.20|) and ()5.7|) the operator P£ in the type A case reads 

P^^(-lfm( q )^z'( q f^, (5.20) 

and thus Vjj- = kerP^ is also the type A monomial space (jf>.l)) . Hence, the solvable sectors 
of the type A Hamiltonians H^- are given by 

V^ = e-<(l,.,...,^- 1 >|_ (g) , (5.21) 

where Wj^- are denned in Eq. ()4.27j) . 

Another important consequence of the fact that the type A systems with arbitrary PDM 
are obtained from the same gauged Hamiltonians (J5.2)) is that we can obtain completely 
the same results as in the constant mass case if we follow the procedure of generating the 
generalized Bender-Dunne polynomial (GBDP) systems {7r^'(£')}^ =0 [3^, 37]. Hence, the 



anti-commutator of the type A A/"-fold supercharges in the PDM case is also proportional 
to the A/th critical GBDP in the type A super-Hamiltonian (for details see Ref. pi) and 
references cited therein): 

{Q M ,Q^}=2 Af ^ ] (U M ). (5.22) 



VI. CLASSIFICATION OF THE MODELS 

Owing to the fact that the form of the gauged Hamiltonians does not depend on the mass 
function m(q), we can classify type A A/"-fold supersymmetric PDM systems completely the 
same way as in the constant mass case In the complex classification scheme based on 
the GL(2,C) invariance, there are five inequivalent models according to different canonical 
forms of A(z) given in the second column of Table [I] 

A new feature due to the position dependence of mass emerges through the change of 
variable. From Eq. ()4.10|) the change of variable is determined by 



±u(q) = ±f d( l V™(q) = J —, ==. (6.1) 



Expressing the variable z in terms of u from this formula, z = f(u), we obtain the change 
of variable as 

z = z{q)=f(u)\ (6.2) 
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The function f(u) in each of the five cases is given in the third column of Table |1J When 
m(q) = 1, it is evident that u(q) = q (up to an additive constant), and thus reproduces 
all the type A models with constant mass in Ref. j^. The form of the change of variable 
(16. 2|) indicates that it is more convenient to express the type A PDM Hamiltonians (|5.13j) in 
terms of f(u). To this end, we must first express derivatives of z with respect to q in terms 
of f(u) and m(q). For instance, the first derivative of z(q) reads 



z'{q) = f'{u)u'{q) = m(g)a/'(u). 



Similarly, we can derive 



z\q)=m{q)f\u) + ^ J f\u), 

2m(q) 2 



z"\q) = m(q)lf'"(u) + lm'(q)f"(u) + ( ^[) _ ^M_)f {u 



\2m(q) 2 4m(g) 



(6.3) 

(6.4) 
(6.5) 



Using Eqs. (jPjj) and (fB~3"jl - (j6~Kjl . we obtain, 

2z"'(q) 3z"(q) 2 



2E'{q) - E(qY 



z'(q) z'(q) 2 

2f"'(u) 3/» 2 



m(q) 



f'{u) f'{uf 



m"(q) 5m'(g) s 



m(q) 4m(q)' 



(6.6) 



The function W(q) and its derivative are also expressed in terms of f(u) and m(q) with the 
aid of Eqs. (JOfijl . IjCTl and (JH3D as 



W(q) 
W'(q) 



/(«) 

m(q)f"(u 



*=/(«) 



m'(q) 



f'(u) 2 2m(g)3/'(u) 



Q(z)-m( g )Q'(z) 



(6.7) 
(6.8) 



«=/(«) 



Substituting Eqs. ()6.6j) - (j6.8j) into Eq. (|5.13jl . we finally have the expression of the type A 
PDM Hamiltonians in terms of f(u) and m(q) as follows: 



d 2 m'(q) d 



2m(q) dq 2 m(q) 2 dq 



V'{u) 



N 2 -l(2f"{u) 3/» 
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_ 7m'(g) 2 AT //>) _ 
+ 8m(g) 2 32m(g) 3 2 Vf («) 2 



Similarly, the gauge potentials Wj^ are expressed as 



/'(« 
- i? 



/'(«) 



(6.9) 



«=/(«) 



--m|m(g)| + — ; 



ln|/'00l± / dw 



WOO) 

/'(«) ' 



Hence, the solvable sectors of the type A Hamiltonians (|5.21|) reads 



V±=m(q)y\u) 



exp =)= / du 



WOO) 

/'00 



(1, /(«),... ./(u)"- 1 ). 



(6.10) 



(6.11) 
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Observing the above formulas and noting the fact that u(q) — > q (up to an additive constant) 
as m(q) — > 1, we can find a simple procedure to obtain a type A PDM quantum system from 
a type A constant-mass model. That is, if we denote the pair of potentials, gauge potentials, 
and solvable sectors of a type A constant-mass model as V^ ± (q), Wj$ (<?), and Vj$ [?], 
respectively, those of the corresponding type A PDM model are given by 

Wfr{q) = -\ln\m(q)\ + W^^q)), Vfr[q] = m(g)*VjJ }± K 9 )]. (6.12b) 

This result is consistent with the one obtained by the point canonical transformation, see 
e.g., Eqs. (2.7) and (2.8) in Ref. @, and Eqs. (10), (13) and (14) in Ref. 0. Advantages of 
the framework of A/"-fold supersymmetry are that it can automatically determine the general 
functional form of the potentials for which we can obtain (a number of) analytic solutions 
and that, in a good situation like the present type A case, we can completely classify all 
the possible (quasi-) solvable models which have a specific type of solutions. In addition, it 
enables us to obtain simultaneously a pair of almost isospectral Hamiltonians. 

In what follows, we show the explicit form of effective potentials and solvable sectors in 
each case classified as in Table HI We note that type A models become not only quasi-solvable 
but also solvable when the parameters in the potentials satisfy 

ct4 = a 3 = 6 2 = 0. (6.13) 

Hence, the models in Cases LTV are solvable when b 2 = 0, while the model in Case V is 
always only quasi-solvable irrespective of the values of the parameters 6j. 

A. Case I: A{z) = 1/2, f(u) = u 
Effective potentials: 

uik) - \ ("Mi? + M,) + M 2 t« + g£ - - R . ( 6 . 14 ) 

Solvable sectors: 



(6.15) 



y± = m(q)' exp (^f^u(q) 3 =F ju(q) 2 T b u(q)^j (l, u(q), u(q) N 1 ). 

B. Case II: A(z) = 2z, f(u) = u 2 
Effective potentials: 

U^{q) = h ^u{qf + b -^u(q) 4 + i(6? + 2b b 2 T QMb 2 )u{qf 

{M-l±b ){M+l±b,) m"(q) _ 7m' (q) 2 Ub x b b x 

8u(q) 2 8m(q) 2 32m(g) 3 + 4 4 ' { ' j 

Solvable sectors: 

= m{q)*u{q)-^ exp (t|«(<?) 4 ? \u{qf^ <1, u(q) 2 , u{qf^). (6.17) 
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C. Case III: A(z) = 2vz 2 , f(u) = e 2 ^ u 
Effective potentials: 

U±( q ) = % e W*>{<i) + b.( bl T 2A/V)e 2 ^) + ^(h ± 2Mu)e^ u ^ + ^ e ~ A ^ u ^ 
A/v ' 8v 4z/ v ' 4z/ v ' 8v 

m"(q) 7m'(q) 2 b\ + 2b 2 b J\f 2 - 1 
+ 8m(g) 2 32m(g)3 + 81/ + 6 1 ' 

Solvable sectors: 

V±=m{q)\ exp (t^ 2 ^ ± ^e" 2 ^® _ * b \ {q) ^ 

D. Case IV: A(z) = 2u{z 2 - 1), f(u) = cosh2^/Fu 
Effective potentials: 

TT ±, \ h l ■ u 2o r i \ , ^iT^) . , 2 ^ , (62 + 6o)(fei±2A^) 

= smh 2 vM<?) + 5- smh V^u(q) + 2 

w Iv 81/ smh ^vu\q) 

(62 + 60 - 61 =F 2(AT - l)i/)(6 2 + 6 - 61 =F 2(AT + 



+ 



^sinh 2 2y / z/w(g) 



7m'(g) 2 ^A/~fe 2 2fe 2 (fe 2 + fe + fe 1 ) + fe 2 A/" 2 - 1 
+ M^"32^F T ^~ + 8^ + _ — (6 ' 20) 

Solvable sectors: 

Vjy =m(g)2(sinh2v / ^'u(g))~^2 _=F ^ ; (tanh v / ^'w(<?)) =F ^ iira exp ^=F^- cosh 2y / ^'M(g)^ 

x <l,cosh2 v ^M(g),...,(cosh2 v ^M(g)) A ''" 1 ). (6.21) 

E. Case V: A(z) = 2z 3 - g 2 z/2 - g 3 /2, f(u) = p{u) 
Effective potentials: 

m"(q) _ 7m'(q) 2 Nh b 2 b, 

8m(q) 2 32m(g) 3 4 4' 1 ' ' 

where e\ = p(wi) (I = 1,2,3) are the values of the Weierstrass function at the half of the 
fundamental periods 2wi and Hf = 3e 2 — g 2 /A. The coupling constants r]f are given by 

•nf = - b 2 ei{b 2 ei - 2b 1 )(2H? - he 2 ) + (b\ + 2b 2 b )e 2 - 2&i&oe, + b 2 

+ (A/" 2 - \){Hf - \8e 2 H 2 + 36ef) t 2Af[(b 2 e t - h){hH 2 - \2e 2 )e l - b Q H 2 } . (6.23) 
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Solvable sectors: 

3 SS-l fyef-HH+bQ 

V£ =m(q)*H\p(u(q)) - ef 4 «T (l, p(u(q)), . . . , pHq))^ 1 ). (6.24) 
i=i 



VII. DISCUSSION AND SUMMARY 

In this paper, we have generalized A/"-fold supersymmetry in ordinary quantum systems 
to those with position-dependent mass. The significant properties such as the equivalence to 
weak quasi-solvability also hold in the PDM case. We have developed the general procedure 
to construct an A^-fold supersymmetric PDM system and applied it to obtain the general 
form of type A A^-fold supersymmetry in PDM quantum systems, it turns out that the 
framework of A/"-fold supersymmetry is quite powerful also in searching (quasi-) solvable 
PDM Hamiltonians as well as ordinary ones. In fact, many of the so far constructed (quasi- 
solvable PDM Hamiltonians in the literature are realized as special cases of type A M- fold 
supersymmetry. In addition, we can simultaneously obtain a pair of almost isospectral PDM 
Hamiltonians in the framework of A/"-fold supersymmetry. 

There are a lot of future issues worthy investigating in this research direction. For exam- 
ple, it would be straightforward to construct other types of A/"-fold supersymmetry, namely, 
type B j^H and type C |3(| in PDM quantum systems. They are characterized by respec- 
tively the type B monomial space 

VV = <M,...,^- 2 ,^>, (7.1) 

and by the type C monomial space 

Vj* = {l,z,...,z">-\z*,z»\...,z»"*- 1 ), (7.2) 

with A G K \ {— A2, —A/2 + 1, • • • ,A/i}. Applying the algorithm in Section |TV| with above 
spaces, we will obtain the general form of type B and type C A/"-fold supersymmetric PDM 
quantum systems. We expect the same relations as Eqs. ()6.12|) since they are universal in 
view of the point canonical transformation. 

Investigation of dynamical properties are interesting and important. In contrast to the 
ordinary constant-mass case, not only the form of potential but also the position dependence 
of mass affect various aspects of PDM systems. In particular, dynamical A/"-fold supersym- 
metry breaking can take place through the nonperturbative effect due to quantum tunneling 
39l 4p| . Hence, it is quite interesting if we can experimentally observe such a phenomenon 
in realistic systems such as semiconductors, quantum dots, and so on. 

Construction of (quasi-)solvable quantum many-body systems with position-dependent 
mass is a challenging problem. Recently, we have developed systematic and powerful meth- 
ods for constructing quasi- solvable differential operators of arbitrary number of variables 
and applied them to construct quasi-solvable many-body Hamiltonians with constant mass 
We expect that the methods would also work well in the PDM case. Up to now, 



we have found four different types of quasi-solvable differential operators of arbitrary number 
of variables, namely, type A j4lj], type C 42||, type A', and type C j43[. In the constant-mass 



case, most of the obtained Hamiltonians are Calogero-Sutherland and Inozemtsev models 
associated with classical root systems. Hence, we will obtain a set of mass-deformed quan- 
tum systems of these types which preserve (quasi-) solvability. In this respect, we anticipate 
a many-body generalization of the relations ()6.12|) . which relates a (quasi-) solvable potential 
in PDM systems with that in constant-mass systems. 



16 



Acknowledgments 



This work was partially supported by the National Science Council of the Republic of 
China under grant No. NSC-93-2112-M-032-009. 



[1] L. Dekar, L. Chetouani, and T. F. Hammann, J. Math. Phys. 39 (1998) 2551. 
[2] L. Dekar, L. Chetouani, and T. F. Hammann, Phys. Rev. A 59 (1999) 107. 
[3] V. Milanovic and Z. Ikonic, J. Phys. A: Math. Gen. 32 (1999) 7001. 

[4] A. R. Plastino, A. Rigo, M. Casas, F. Garcias, and A. Plastino, Phys. Rev. A 60 (1999) 4318. 
[5] A. de Souza Dutra and C. A. S. Almeida, Phys. Lett. A 275 (2000) 25. quant-ph/0306065. 
[6] A. D. Alhaidari, Phys. Rev. A 66 (2002) 042116. quant-ph/0207061. 

[7] B. Gonol, B. Gonol, D. Tutcu, and O. Ozer, Mod. Phys. Lett. A 17 (2002) 2057. quant- 
ph/0211112. 

[8] B. Gonol, O. Ozer, B. Gonol, and F. Uzgiin, Mod. Phys. Lett. A 17 (2002) 2453. quant- 
ph/0211113. 

[9] R. Koq, M. Koca, and E. Korciik, J. Phys. A: Math. Gen. 35 (2002) L527. quant-ph/0410108. 
[10] B. Roy and P. Roy, J. Phys. A: Math. Gen. 35 (2002) 3961. 

[11] A. de Souza Dutra, M. Hott, and C. A. S. Almeida, Europhys. Lett. 62 (2003) 8. hep- 
th/0306078. 

[12] R. Koc and M. Koca, J. Phys. A: Math. Gen. 36 (2003) 8105. quant-ph/0410127. 

[13] R. Koc and H. Tiitunciiler, Ann. Phys. (Leipzig) 12 (2003) 684. quant-ph/0410088. 

[14] K. A. Samani and F. Loran. Shape invariant potentials for effective mass Schrodinger equation. 

Preprint, quant-ph/0302191. 
[15] J. Yu, S.-H. Dong, and G.-H. Sun, Phys. Lett. A 322 (2004) 290. 
[16] J. Yu and S.-H. Dong, Phys. Lett. A 325 (2004) 194. 
[17] G. Chen, Phys. Lett. A 329 (2004) 22. 
[18] G. Chen and Z. Chen, Phys. Lett. A 331 (2004) 312. 

[19] C. Quesne and V. M. Tkachuk, J. Phys. A: Math. Gen. 37 (2004) 4267. math-ph/0403047. 
[20] B. Bagchi, P. Gorain, C. Quesne, and R. Roychoudhury, Mod. Phys. Lett. A 19 (2004) 2765. 
quant-ph/0405193. 

[21] B. G6n61 and M. Kocak, Chin. Phys. Lett. 22 (2005) 2742. quant-ph/0412161. 

[22] B. Bagchi, A. Banerjee, C. Quesne, and V. M. Tkachuk, J. Phys. A: Math. Gen. 38 (2005) 

2929. quant-ph/0412016. 
[23] S.-H. Dong and M. Lozada-Cassou, Phys. Lett. A 337 (2005) 313. 
[24] B. Roy and P. Roy, Phys. Lett. A 340 (2005) 70. 

[25] B. Bagchi, P. Gorain, C. Quesne, and R. Roychoudhury, Europhys. Lett. 72 (2005) 155. quant- 
ph/0505171. 

[26] C. Quesne. First-order intertwining operators and position-dependent mass Schrodinger equa- 
tions in d dimensions. Preprint, quant-ph/0508216. 
[27] O. von Roos, Phys. Rev. B 27 (1983) 7547. 

[28] A. A. Andrianov, M. V. Ioffe, and V. P. Spiridonov, Phys. Lett. A 174 (1993) 273. hep- 
th/9303005. 

[29] H. Aoyama, M. Sato, and T. Tanaka, Nucl. Phys. B 619 (2001) 105. quant-ph/0106037. 
[30] A. Gonzalez-Lopez and T. Tanaka, J. Phys. A: Math. Gen. 38 (2005) 5133. hep-th/0405079. 



17 



A. V. Turbiner and A. G. Ushveridze, Phys. Lett. A 126 (1987) 181. 

A. G. Ushveridze, Quasi-exactly Solvable Models in Quantum Mechanics (IOP Publishing, 
Bristol, 1994). 

H. Aoyama, M. Sato, and T. Tanaka, Phys. Lett. B 503 (2001) 423. quant-ph/0012065. 
T. Tanaka, Nucl. Phys. B 662 (2003) 413. hep-th/0212276. 
A. V. Turbiner, Contemp. Math. 160 (1994) 263. 

A. A. Andrianov and A. V. Sokolov, Nucl. Phys. B 660 (2003) 25. hep-th/0301062. 

C. M. Bender and G. V. Dunne, J. Math. Phys. 37 (1996) 6. hep-th/9511138. 

A. Gonzalez-Lopez and T. Tanaka, Phys. Lett. B 586 (2004) 117. hep-th/0307094. 

H. Aoyama, H. Kikuchi, I. Okouchi, M. Sato, and S. Wada, Nucl. Phys. B 553 (1999) 644. 

hep-th/9808034. 

M. Sato and T. Tanaka, J. Math. Phys. 43 (2002) 3484. hep-th/0109179. 

T. Tanaka, Ann. Phys. 309 (2004) 239. Erratum-ibid. 320 (2005) 257, hep-th/0306174. 

T. Tanaka, Ann. Phys. 316 (2005) 187. hep-th/0407275. 

T. Tanaka, Ann. Phys. 320 (2005) 199. hep-th/0502019. 



18 



